Symmetry, Integrability and Geometry: Methods and Applications 



Vol. 2 (2006), Paper 073, 11 pages 



Quasi-Exactly Solvable AT-Body Spin Hamiltonians 
with Short-Range Interaction PotentialsQ 

A. ENCISO, F. FINKEL, A. GONZALEZ-LOPEZ* and M.A. RODRIGUEZ 
Depto. Fisica Teorica II, Universidad Complutense, 28040 Madrid, Spain 



E-mail: artemio@fis.ucm.es 



Received September 15, 2006, in final form October 23, 2006; Published online November 03, 2006 



Original article is available at http://www.emis.de/journals/SIGMA/2006/Paper073/ 



Abstract. We review some recent results on quasi-exactly solvable spin models presenting 
near-neighbors interactions. These systems can be understood as cyclic generalizations of 
the usual Calogero-Sutherland models. A nontrivial modification of the exchange operator 
formalism is used to obtain several infinite families of eigenfunctions of these models in 
closed form. 
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1 Introduction 

In the early 1970s, F. Calogero [6] and B. Sutherland [37] introduced the quantum integrable 
systems that nowadays bear their names. Apart from their intrinsic mathematical interest \25\ 
EH [3], Calogero-Sutherland (CS) models have played a central role in Physics due to their 
relevant applications to as diverse topics as soliton theory [231 El], quantum field and string 
theory [181 [9] , quantum Hall effect [2] , fractional statistics [28] and random matrix theory [36] . 
The first satisfactory explanation of the integrability of these models was given by Olshanetsky 
and Perelomov [26] , who connected these models with the radial Laplacian of symmetric spaces 
associated to the root system A]sf. This unified view enabled them to introduce important 
generalizations, including different root systems and elliptic potentials. 

During the last decade CS models have been extended to the case of particles with internal 
degrees of freedom, which we shall henceforth call spin. There are two different approaches to 
spin CS models, namely the supersymmetric [5] and Dunkl operator \10\ |29] formalisms. These 
methods have allowed to solve totally or partially several rational, trigonometric and elliptic 
spin models, both in their Aj^ and BCn versions [8| I42[ [TBI I17j . The interest in spin CS models 
also increased as a consequence of their direct connection with the Haldane-Shastry (HS) spin 
chain |19[ , which was laid bare by Polychronakos through the so called "freezing trick" [30] . 
This technique was also used in the construction of solvable spin chains associated to different 
potentials and root systems [3ll [l2] . 

Auberson, Jain and Khare |22[ [T] introduced partially solvable versions of the CS models 
in which each particle only interacts with its nearest and next-to-nearest neighbors. Similar 
scalar models were also studied by Ezung, Gurappa, Khare and Panigrahi p2|- There are two 
reasons that make these kind of systems very promising from a physical point of view. First, 
some of them are related to the short-range Dyson model in random matrix theory [3j. Second, 

*This paper is a contribution to the Proceedings of the Workshop on Geometric Aspects of Integ- 
rable Systems (July 17-19, 2006, University of Coimbra, Portugal). The full collection is available at 
|http: / /www. emis.de/journals / SIGMA/Coimbra2006.html , 
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the HS chains associated to these models would occupy an interesting intermediate position 
between the Heisenberg chain (short-range, position-independent interactions) and the usual 
HS chains (long-range, position-dependent interactions). A first step towards the construction 
of these chains was Deguchi and Ghosh's definition of spin 1/2 versions of the Jain and Khare 
Hamiltonians [7] using the supersymmetric formalism. Unfortunately, all these authors solely 
managed to construct a few exact solutions, and all of them with trivial spin dependence, and 
the procedures developed to obtain exact solutions are by no means systematic. 

In |13^ ll4j we introduced three new families of spin near-neighbors models and used a nontriv- 
ial modification of the Dunkl operator method to obtain a wide range of fully explicit solutions. 
In this article we aim to review the main ideas underlying our constructions and explain their 
connection with the usual CS models. In Section [2] we define the Hamiltonians that we shall deal 
with and state our main result, which is a description of the algebraic states of the models. In 
Section [3] we review the philosophy underlying the calculation of exact solutions of Schrodinger's 
equation by algebraic methods. In Section H] we sketch the main logical steps that the construc- 
tion of the invariant flags rests upon. We finish the paper by showing in Section [5] how the 
actual computation of the algebraic eigenfunctions can be carried out. For the sake of brevity, 
complete proofs are not given in this review. 



2 Main result 

Let S be the Hilbert space of internal degrees of freedom of particles of spin M G iN. Let 
us fix a basis 

B={\su...,sn) : Si £ {-M,-M + 1,...,M}} 

of S and define the spin exchange operators Sij as 

Sij I Si , • • • , Sj , . . . , Sj , . . . , Stv) — I Si , . . . , , . . . , Sj , . . . , s^) . 

These operators can be expressed in terms of the (normalized) fundamental generators of 
SU(2M + 1) as Sij = (2M + 1)"^ + J^J^, the index sum ranging from 1 to 4M(M + 1). 
The Hamiltonians of the models we shall be concerned with are given by (cf. |14j ) 

He = -Y,dl + V,, 6 = 0,1,2, (1) 

i 

where 

Vo = + Y: 7 ^4 T + E ^^^(« - (2a) 

^ [Xi - Xi-l){Xi - Xi+i) ^[Xi-Xi+iy 

h{h - 1) >^ Sa^xf 



xf + x'^ 



r\2 (^-w), (2b) 

V2 = 2o^^cot(xi - Xj.i) cot(xi - Xi+i) + 2a^csc^(xj - Xi+i){a - Si^i+i), (2c) 



with = xf and a,b > 1/2. Here and in what follows, all sums and products run from 1 to A^ 
unless otherwise stated, with the identifications xq = xjy and xjsi^i = xi. There is a hyperbolic 
potential analogous to (pcj) that is recovered substituting Xi by — Xj and V2 by — V2. The scalar 
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reductions H^'^ of these models are obtained from the above Hamiltonians by the substitution 
Si^i+i — > 1- 

A few remarks on the configuration spaces of these models are now in order. In all three 
models the potential diverges as (xj— on the hyperplanes Xj = Xj+i, so that the particles i 
and i + 1 cannot overtake one another. Since we are interested in models with nearest and next- 
to-nearest neighbors interactions, we shall henceforth assume that xi < ■ ■ ■ < rrjv. For the 
second potential ()2bp we shall take in addition xi > 0, due to the double pole at Xj = 0. 

A first observation concerning the eigenstates of the models ([T]) is that if ^ is an eigenfunction 
of H^'^ with energy E, then the factorized state ^ = tpls) is an eigenfunction of ff^ with the same 
energy for any spin state \s) symmetric under permutation of particles. A second observation 
is that H2 commutes with the total momentum P = — i^li^i^i- Hence the movement of the 
center of mass decouples, and only the eigenfunctions of H2 with zero total momentum need to 
be considered. 

The next theorem summarizes the main results presented in |13l I14j . which yield a fully 
explicit description of several families of algebraic eigenfunctions. It is not difficult to realize, 
however, that these eigenfunctions do not exhaust the whole spectrum of the models. We 
shall denote by A the projection operator on states totally symmetric under the simultaneous 
permutation of both the spatial and spin coordinates. We define the spin vectors \si), \sij) as 



A(xi|s)) = ^Xi|Sj), A(xiX2|s)) = y^^XjX j\Si 

i i<j 



We shall consider the subspace S' C S of spin vectors \s) such that Ylil^i) symmetric. 
A thorough characterization of this space is given in 



Theorem 1. Let I, m he nonnegative integers and let \ s) be an arbitrary spin vector, and denote 
by X the center of mass coordinate Xj. Then the following statements hold: 

1. Let a = N{a + i) - |, /3 = /3(m) = 1 - m - N{a + t = ^ - 1, and /i/^ = 
Xj — Xj-i-il'*. The Hamiltonian Hq possesses the following families of spin 
eigenfunctions with eigenvalue Ei^ = Eq + 2uj{2l + m), where Eq = Nuj{2a + 1) is the ground 
state energy: 



fl=l^lmPl^^"m\s), m>0. 



plZtl:^\t)K{x,\s)) - xPlZtlf\t)K\s) 



T (2) 



>{"+2,/3) 



f (t)(A(x2|s))-2xA(xi|s))) 



+ x^ 



(3) 



'■^l'\t)Y,4 + x\m{t) A|S), 
i 

^ i-x2p(^+Y^(t)A(x?|5)) 



2([ 



m— 3 ] _| 

2 J 2^ 



i 



m > 1, 



m > 2, 



m > 3, 



m > 4. 
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When \s) G S' there is an additional family of eigenf unctions given by 
Pll^lf\t){H^iX2\s)) - 2xA(xi|s))) 



T(2) 



T/ie functions (fm, (pm o,nd Xm o,re polynomials given explicitly by 



m>2. 



_ m + 2a + 2 (o+2,/3-2) _ p(Q+3,/3-i) _ 4a + 7 (a+2,/3-1) 1 p(a+3,/3) 

O 2 



m — 1 



m — 1 2 



{«+4,/3-l) _ 2p(a+3,/3-l) _ m + 2a + 3 p(a+2,/3-l) 



^ rr 



{rn — l){in — 3) 2 ^ 



1 p{a+4,/3) m + 2a - 1 (q+3,/3) 

62 m, — 6 2^ 

3m + 2a (a+2,/3-1) 2m - 7 p(a+3,/3-i) _ p{a+4,/3-i) 



(m — l)(m — 3) 2 



m _ -1 
n J- 



m — 3 



m _ 1 



(m — l)(m — 3) 2 
/or ewen m, and 



m + 2a + 2 (a+2,/3) _ m, + 2a (a+3,/3) 1 p{a+4,/3) 



m — 3 



,^ _ 9 p(a+2,/3-l) _ p(a+3,/3-l) 1 p(a+3,/3) + + 2 (a+1,/3) 

3 — 2~ m(m — zj — 2~ 



2 2 

m + 2a + 2 



-P 



m — 3 5 
2 



m — 2 

, _ p(a+4,/3-l) _ 2m - 5 {«+3,/3) _ 1 p(o+4,/3) m + 2a - 1 (a+3,/3) 

O 2 



m — 2 



m — 2 



m — 5 1 
2 



2m - 3 (a+2,/3-1) 2(m - 3) (a+3,/3-1) _ (a+4,/3-1) 



m{m, — 2) 



2- m — 3 
— 2~ 



m — 3 
2 



m + 2a + 1 p(Q+2,/3) _ m+_2a («+3,/3) 1 (a+4,/3) 

m — 2 3 — 2~ 



m(m — 2) 
/or oc/d m. 

2. T/ie Hamiltonian Hi possesses the following families of spin eigenf unctions with eigenvalue 
Ek = Eq + Aku), where Eq = Nuj{4:a + 26+1) is the ground state energy: 

,(0) _ ..ra~l/, .J2\ 



(1) 
k 

,(2) 



^lLl~'{uJr')^\s), k>0, 
^iLltli^r^) [Nk{xl\s)) - r^K\s)] , k>l, 



loii/i a = A^(2a + 6 + i), /? = iV(4a + 6 + |) and fi = e"'^'"' Hi l^^i - ^f+iT^;- . 

3. T/ie Hamiltonian H2 possesses the following spin eigenf unctions with zero momentum 



*3 = ^ 



2a 



2a + 



_$(0) + ^cos(2(x,-Xj))|sj 



= /i^sin(2(xi - Xj))|sj), 



where /u = Hi^i'^'' ~ Their energies are respectively given by 

1 



Eo, Eio = EQ + A[2a + l 



^3,4 = ^o + 8(2a + 1), 



where Eq = 2Na^ is the ground state energy. 
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3 Calogero models and Dunkl operators 

Let us consider a self-adjoint operator H acting on a given Hilbert space 7i. The idea underlying 
the construction of algebraic eigenfunctions of H (see, e.g., [lUj) is that the explicit knowledge 
of a finite-dimensional subspace Tii C TC which is invariant under the operator H allow ones to 
compute dim 7^1 eigenfunctions and eigenvalues of H by algebraic methods, i.e., diagonalizing 
the matrix -ffl'Hi- In this case the operator is said to be quasi- exactly solvable (QES). In fact, 
the Hamiltonians of the models (j2ap and ()2bp possess an infinite flag 7i\ C 7i2 C • • • of known 
finite-dimensional invariant subspaces, which yield an arbitrary large number of eigenvalues and 
eigenfunctions. Although such models are sometimes termed "exactly solvable" |39j . we will not 
use this terminology since the algebraic eigenfunctions of the models ()2ap and ()2bp are not an 
orthonormal basis of the Hilbert space. 

A particularly convenient method of carrying out this program |381 [35] is through a Lie 
algebra g of first-order differential operators J" = £,°'^{z)dzi + ^°(^) = 1) • • • j ^) acting on 
a finite-dimensional module M C C°°(M), with M a domain in M.^ . Let us assume that there 

^ r r 

exists a second-order differential operator H = ^ CahJ°'J^+^ CaJ"" + cq that is equivalent, 

a,b=l a=l 

up to gauge transformation H ^ ^Hfi~^ [fi G C°°(M, M"*")) and global change of variables z G 
to a Schrodinger operator H = — + V , standing for the Laplace-Beltrami 
operator in (M^,^). If ^iM\z^x C L"^ {M.^ , ^dx) , then one can obtain dimA^ eigenfunctions 
of H by algebraic methods. The Schrodinger operators amenable to this treatment are termed 
Lie-algebraic. It should be observed that the operators J" are not symmetries of H: q is an 
algebra of hidden symmetries of H. 

Dunkl operators [20] were originally introduced to study spherical harmonics associated to 
measures invariant under a Coxeter group [10]. Actually, let v,z G and define the reflection 
a^z = z — 2\v\~'^{z • v)v. Let 2B be a Coxeter group, which can be assumed to be the Weyl 
group of a (possibly nonreduced) root system R, and let Va = vl^Si (a = 1, . . . , r) be a basis 
of positive roots. Define an action of 2IJ on C[z] as Kaf = f o (^va- Dunkl operators were 
originally defined as 

J^ = ^.^+Y,^{l-Ka), (3) 

where the real parameters are chosen so that they are constant on each orbit of 2IJ. The 
operators Jj can be understood as deformations of the partial derivatives dz^ that commute 
with the deformed Laplacian, i.e., Yli Ji-^j — Jj Si Jf- It can be verified that { Jj, Ka} span a 
degenerate Hecke algebra j24j . Currently, the definition of Dunkl operators has been generalized 
to mean a set { Jj}^^ C EndC[2;] of first order differential operators which leave invariant finite- 
dimensional polynomial subspaces and such that {Ji,Ka} span a Hecke algebra. 

In their original form, Dunkl operators are directly connected with the rational Calogero 
model of type R [H]. In the simplest case {R = ^iv-i, 2U = Sn) and writing x instead of z, 
the Dunkl operators read 

Ji = d,^+g^y2^—{l-K,j), 

where 

{Kijf){xi, ...,Xi,...,Xj,.. .,xn) = /(xi, ...,Xj,...,Xi,.. .,xn) (4) 
denotes the reflection operator associated to the root ej — Cj. Actually, let us consider the 
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Calogero ground state function /x(a;) = e~^^ Tli<j ~ ^^'^ the auxiliary operator 

.f = ^x,d,^. (5) 

i 

A straightforward calculation shows that one can use the ground state function fx to gauge 
transform the deformed Laplacian so that 

^^Hcfi-' = -^dl+aY^ ^ (a - Kij) + u;r^ (6) 
with = a(a — 1) and 

i 

yields the Calogero model of Aj^-i type when acting on symmetric functions ip E AL^(M^). 
Since Jj, preserve the space of polynomials 

= {/gC[^] :deg/<7^} (7) 

for any n = 0, 1, ... , the Dunkl operators provide a very convenient fashion of exploring the 
solvability properties of Calogero-Sutherland models. It should be remarked, however, that for 
an arbitrary Coxeter group 23J the Dunkl operators ([3]) do not form a Lie algebra; nevertheless, 
this technique captures most of the relevant features of the Lie-algebraic method. 

As we shall now outline, Dunkl operators can also be used to introduce internal degrees of 
freedom in the picture without breaking the solvability properties of the models. Given a scalar 
differential-difference operator D linear in Kij, let us denote by D* the differential operator 
acting on C°° S obtained from D by the replacement Kij — > Sij. It is clear that the actions 
of D and D* coincide on the (bosonic) Hilbert space 

W = A(l2(M^) ® S). (8) 

Therefore the scalar operator ([6]) coincides with that of the spin Calogero model 

when acting on symmetric states ^ £ 7i. 



4 Invariant subspaces 

The proof of the main theorem rests on the construction of appropriate invariant flags for the 
Hamiltonians ([1]) using a modification of the Dunkl operator formalism. This modification turns 
out to be rather nontrivial, ultimately due to the fact the cyclic group is not of Coxeter type. 
As the first step, we consider the second-order differential-difference operators given by 

Zi — Zi-^l '—^ \Zi — Zi-^-l) 

where di = dzi, zn+i = zi^ and 



'9o{x,y) = l, 'di{x,y) = -{x + y), §2{x,y)=xy. 
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Each Hamiltonian is related to a linear combination 

H, = cT, + c_ J- + Co J° + Eo (10) 

of its corresponding operator and the auxiliary first-order differential operators J~ = ^ - di 
and = Zidi via a change of variables, a gauge transformation, and the star mapping 
defined in the previous section, that is. 

The constants c, c_, cq, Eq, the gauge factor /x, and the change of variables C for each model 
are listed in Table [TJ Hence the construction of the models ([1]) is analogous to that of the usual 
Calogero-Sutherland models, with the operators ([9]) being a cyclic analog of the sum of the 
squares of the Dunkl operators. 



Table 1. Parameters, gauge factor and change of variable in equations (jlOp and (fTTj) . 





e = 


e = 1 


e = 2 


c 


-1 


-4 


4 


c_ 





-2(26 + 1) 





Co 


2uj 


4u; 


4(1 - 2a) 


Eq 


Nuj{2a + 1) 


Afa;(4a + 26 + 1) 


2iVa2 


Kx) 


e^'t''' n - ^i+ll" 

i 


— — TT 1 2 2 If* h 

e 2- -^ml < 

i 


Hsin" \xi - Xj+il 

i 


C(x) 


X 


2 


g±2ix 



It is obvious that the operators ()10p preserve the polynomial space ([7]) for any n G N, so 
one may be tempted to believe that the usual arguments for spin CS models should yield an 
invariant flag for the Hamiltonians ([T]). Nevertheless, this is not the case. In fact, the standard 
construction |16^ [T7] is based on the fact that the actions of Hq and fiHcfi~^ on symmetric 
states coincide and these operators commute with the Sn symmetrizer A. Unfortunately, 
(or H^) do not commute with A and cyclic symmetry does not suffice to exchange and 
IJ,H^jj,~^, so this procedure does not grant the existence of any nontrivial invariant subspaces, 
not even of direct products M (S> AT, {M C C°°(M^)). 

We shall now review the actual construction of the invariant spaces, which is considerably 
more involved. We shall not provide complete proofs, but merely a sketch of the main logical 
steps the construction rests upon. 

Classical results on the theory of invariants |41j are responsible for the success of studying 
CS models through symmetric polynomials |27^ I33j. We shall extend this approach to deal with 
the Hamiltonians ([T]). Let us first introduce two bases {cTk} and {rk} of the space of symmetric 
polynomials in z: 

(Tk = J2^i^ -Tfc = J2 ^n-'-^ifc! k = l,...,N. 

i ii<-<ik 

The operators consist of three summands which are of second, first and zeroth order in 
the derivatives. Let us denote each summand by L^, 2aX^ and —2aA^ respectively. It is not 
difficult to realize that {L^} span a Lie algebra isomorphic to 5[(2), as in the CS case. It turns 
out that the first-order differential operators leave invariant a flag of symmetric polynomials, 
as stated in the following easy lemma. 
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Lemma 1. For each n = 0, 1, . . . , the operator leaves invariant the linear space , where 
Remark 1. It should be noted that these flags cannot be trivially enlarged, since, e.g., 

i 

-XiU^ = 2(72 + ^ ZiZi+i, XiTN-l = TN ^(ZjZi+i)""'^, 



ix2(T2 = 2CT2 + ^ 2:^2:^+1, X2TN-2 = Ntn-2 " Tat ^( 



ZiZiJ^l) ^ 



2 

i I 
are not symmetric polynomials. 

In the next proposition we characterize subspaces of the flags described in Lemma [1] that 
are preserved by the whole operator T^. If / G C[cri, cr2, cs, 'T'tv-i, tat], we adopt the convenient 
notation 

^ fa^J, A; = 1,2, 3, 
\drj. k = N-l,N. 

Proposition 1. For each n = 0, 1, ... , the operator leaves invariant the linear space 5", 
where 

5o" = {/ G I /33 = 0}, 

S'l = {/ G X^ I /22 = fNN = 0}, 

S2 = {f ^ '^2 I /ll = fN~l,N~l = 0}. 

Proposition [1] implies that each operator preserves product symmetric subspaces 5" (8) AS 
spanned by factorized states. The main result on invariant subspaces shows that in fact the latter 
operator leaves invariant a richer flag of nontrivial finite-dimensional subspaces of A('P'" ® S). 

Theorem 2. Let 

V = (/(0"l>0'2,0-3)A|s), 5(0-1, 0-2, Cr3)A(2i|s)),/l(cri,CT2)A(Zi|s)), 

h{ai,a2)h{ziZ2\s')) \ /as = 533 = O), 

'^l = (/(0-l,'72,r7v)A|s),g(cri,TAr)A(zi|s)) | /22 = JnN = 9NN = O), 

'^2 = (/(o-i,r7v„i,rAr)A|s),g(TAr_i,rAr)A(zi|s)),TArg(o-i,rAr)A(zf^|s)) 

I /ll = fN-i,N-i = 9n-i,n-i = qii = 0), 

where \s) G S, \s') G S', deg/ < n, degg < n — 1, deg/i < n— 2, degh < n — 2, degq < n—N+1, 
and deg is the total degree in z. Then is invariant under for all n = 0,1, ... . 

From this theorem one easily obtains the following corollary, which is crucial for the compu- 
tation of the algebraic eigenfunctions of the models ([T]). 

Corollary 1. For each e = 0,1,2, the gauge Hamiltonian leaves invariant the space 
defined by 

/^o = ^ ' ^1 = ^1 |/jv=<;^=0' ri2 = h- (12) 
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5 Spectrum and eigenf unctions 

As happens with the usual CS models, the algebraic eigenvalues of the Hamiltonians (pQ) can 
be easily obtained by choosing a basis of the invariant spaces fiTC^ in which the action of H^^ is 
triangular. In all three cases, the algebraic eigenvalue £"0 is the ground state energy, since the 
corresponding eigenfunctions do not vanish in the configuration space C^- 

Now we shall outline how the algebraic eigenfunctions in Theorem [T] were calculated. The 
easiest case is e = 2, since states related by a multiplicative factor only differ by the movement 
of the center of mass, which is conserved. Hence the invariant subspaces Ti.2 solely allow one 
to compute five eigenfunctions of zero total momentum, which are the ones listed in the main 
theorem. 

In the case e = 1 the eigenvalue equation reads 
(^1 - Eq- 4wn)$ = 0, 

where $ is a symmetric vector-valued polynomial in Tii of degree n. Setting t = ujai and 
^ = [p{t) + a2q{t)]A\s) + g{t)A{zi\s)), this equation can be written as 



CtTig = Ctl-^q = 0, C^-'p = -—g - j^tq, (13) 
with the Laguerre operator defined as 

{ct/m = tf"{t) + (A + 1 - t)f{t) + uf{t). 

With some effort one can obtain all the solutions of these equations, as shown in the following 
proposition. 



Proposition 2. The polynomial solutions of the system of ODE's (jl3p are spanned by 

$f =Lri(i)Ak), n>0, 

$W =L^l|(t)[iV^A(zi|s))-tA(zi|.))], n> 1, 

^^n^ = + '^^2 - (3t^]Ms), n>2. 

These solutions are characterized by the conditions q = g = 0, q = and g = respectively 
and correspond to the algebraic eigenfunctions of Hi presented in Theorem [TJ 

The case e = is similar, but the computations become more involved due to the rich 
structure of the invariant flag. Writing 

$ = (p + (T3g)A|s) + (n + a3v)A{zi\s)) + hA{zl\s)) + hA{xiX2\s)), 

where deg<I> = k and /i = if |s) S', the equation (Hq — Eq — 2Lok)^ = reduces to the 
systems of PDE's 

[Lo - 2uj{k -2)]h- 8h2 = 0, (14a) 

[Lo - 2uj{k -2)]h- 8h2 = 6v, (14b) 

[Lo - 2uj{k -l)]u- 4u2 = 4/ii + 4hi + 60-2^1 + 6(2a + l)aiv, (14c) 

[Lo -2uj{k- 4)] V - 16v2 = 0, (14d) 

4a - 

(Lo - 2Lok)p = 2ui + 2(2a + l)h - — — -h + 60-2^1 + 6(2a + l)aiq, (14e) 

[Lo - 2uj{k - 3)]q - 12q2 = 2vi, (14f) 
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with 

Lo = - [Ndl^ + Acjid^^d^^ + Aa2dl^ + 2(2a + l)Nd,;) + 2uj{aid^^ + 2a2a^J. 

One can check by inspection that the system (|14p possesses six families of polynomial solutions, 
as collected in Table [21 These correspond to the six families of algebraic eigenfunctions listed in 
Theorem [TJ 

Table 2. The six types of polynomial solutions of the system and their corresponding eigenfunctions. 







Conditions 


Corresponding eigenfunction 


q = 


u = 


v = h = h = 0, p/0 


*!°; 


u = 


V = 


h = h = 0, 




q = 


V = 


h=h = 0, Uy^O 




q = 


V = 


h = 0, hy^O 




q = 


V = 


h = 0, h^O 




h = 


0, 
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